Nature of the Perpendicular-to-Parallel Spin Reorientation in a Mn-doped GaAs 
Quantum Well: Canting or Phase Separation? 
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It is well known that the magnetic anisotropy in a compressively-strained Mn-doped GaAs film 
changes from perpendicular to parallel with increasing hole concentration p. We study this reorien- 
tation transition at T = in a quantum well with Mn impurities confined to a single plane. With 
increasing p, the angle 8 that minimizes the energy E increases continuously from (perpendicular 
anisotropy) to n/2 (parallel anisotropy) within some range of p. The shape of _E mm (p) suggests 
that the quantum well becomes phase separated with regions containing low hole concentrations 
and perpendicular moments interspersed with other regions containing high hole concentrations and 
parallel moments. However, due to the Coulomb energy cost associated with phase separation, the 
true magnetic state in the transition region is canted with < 8 < n/2. 



During the past decade, the ferromagnetic transition 
temperature of Mn-doped GaAs films has been quickly 
approaching room temperature [l[ . Theoretical studies of 
Mn-doped GaAs based on the Kohn-Luttinger (KL) 
model have been quite successful at modeling and pre- 
dicting much of the behavior found experimentally. In 
agreement with theory Q , experiments show that at 
low temperatures, the magnetic anisotropy of films under 
compressive strain is perpendicular or out of the plane 
when the hole concentration p is small, transforming to 
parallel or in the plane as p increases [5J . Although Mn- 
doped GaAs quantum wells have also been studied the- 
oretically for several years @, 0], little is known about 
the nature of the spin reorientation in a quantum well. 
We show that the spin reorientation in a quantum well 
happens in three stages: for small hole concentrations, 
the angle 9 of the magnetization with respect to the film 
normal is so that the magnetization lies perpendicular 
to the plane; for high hole concentrations, = 7r/2 so 
that the magnetization lies in the plane of the quantum 
well. In between, the moments are either canted with 
< 8 < ir/2 or phase-separated with regions contain- 
ing low hole concentrations and perpendicular moments 
interspersed with regions containing high hole concentra- 
tions and parallel moments. 

The magnetic anisotropy of a quantum well sensitively 
depends on the spin-orbit coupling. In pure GaAs, the 
spin-orbit coupling plays two roles. First, it lowers the 
energy of the j = 1/2 band compared to the j = 3/2 band 
at the T point by about 320 meV. Since the lower j = 
1/2 band is rarely occupied by any holes, it is commonly 
ignored. Second, spin-orbit coupling changes the energies 
of the j = 3/2 holes so that heavy (m/i = 0.5m) and light 
(mi = 0.07m) holes carry angular momentum j-k = ±3/2 
and ±1/2, respectively, along their momentum direction 
k. In the absence of elastic strain, the energies of the 
light and heavy hole bands in bulk GaAs are degenerate 
at the r point. 



In a quantum well bounded by z = ±L/2, the square 
of the z component of the momentum is quantized. For 
the two lowest wavefunctions ipi(z) = y2/L cos(7r z/L) 
and ip2(z) = yJljL sin(27rz / L) of the quantum well, 



{n\kj\m) — (nn/LyS n 



Due to the difference between 



the light and heavy hole masses, the confinement of the 
holes in a quantum well breaks the degeneracy of the 
j = 3/2 bands at the T point with k^ = 0. To simplify 
the following discussion, we shall discuss hole rather than 
electron bands so that hole energies increase quadrati- 
cally like k\ for small k± . Including the effects of lattice 
strain, the energy gap between the j z — ±3/2 and ±1/2 
sub-bands of ip n (z) at the T point is 
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where Q e = e zz - (e xx + €y y )/2 and b d w -1.6 eV is 
the deformation potential [4j. Hence, compressive strain 
(Q f > 0) plays qualitatively the same role as carrier 
confinement within the quantum well. As the quantum 
well becomes narrower, the dominant contribution to this 
splitting comes from the confinement of the holes and 
not from strain, which is typically less than 0.5%. For a 
quantum well with 20 layers or fewer, the contribution of 
strain to A„ can be safely neglected. On the other hand, 
GaAs films with A larger than about 50 cannot be treated 
as quantum wells because too many wavefunctions ip n (z) 
would be required. In such films, the main contribution 
to the band splitting comes from strain rather than 
the confinement of the holes. Regardless of their origin, 
the energy splittings A„ produce a gap in the spin-wave 
spectrum that allows a two-dimensional layer of Mn spins 
to order ferromagnetically Q. 

The carriers in the lower, j z = ±3/2 and upper, 
j z = ±1/2 sub-bands have bandmasses m a = 0.089m and 
mb = 0.197m, respectively, where 4/m a = 3/m; + 1/rrih 
and 4/mb = 1/m; + 3/m/i 0, S|- The coupling of the 
j = 3/2 holes with the Mn spins is included by treating 
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FIG. 1: The energies of the bottom of the carrier bands with 
kj_ = versus J c when only ipi (z) is considered. The chemical 
potential increases with the hole filling. 



the 5 = 5/2 Mn spins classically and by assuming that 
Mn impurities with concentration c are restricted to the 
z = plane. We shall denote p as the number of holes for 
each of the N cations in the central plane. For small p, 
the holes occupy a small portion of the kj^ Brillouin zone 
centered around kj^ = 0. Since the holes then interact 
with many different Mn moments, the precise locations 
and structure of the Mn impurities are not important and 
their interactions with the holes may be treated within 
a mean-field approximation. The exchange coupling of 
the Mn spins Sj = Snij with the holes is then given by 
V = — 2 J c J2i m i • ji, where ji are the hole spins and the 
sum is over all Mn sites. Comparing V with the potential 
used in Refs.[3[ and [4|, we obtain the exchange coupling 
J c = (S/4j)(3N 2c/A, where /3N w 1.2 eV is estimated 
from photoemission measurements [3] and A = L/b is the 
number of layers in the quantum well. Here, b w AA is 
the Ga lattice constant in the z — plane. It follows that 
J c « 1 c/A eV. 

While our results are supported by numerical calcula- 
tions that include both wavefunctions tpi(z) and ip%(z), 
a qualitative understanding of the magnetic anisoptropy 
can be obtained from a simplified model that considers 
only ipi( z )- Due to the typically small Mn concentra- 
tions, the demagnetization fields can be neglected com- 
pared to the anisotropy introduced by the electronic band 
structure. If the Mn moments are tilted an angle 8 away 
from the z axis, then to linear order in J c , the two lower 
sub-bands are split by ±3J C (9 — 0) or (9 = 7r/2), as 
shown in Fig.l for A = 10. By contrast, the two up- 
per bands are split by ±J C (9 = 0) or ±2J C (9 = tt/2). 
When only the lowest sub-band is populated by holes, 
the energy difference (E(9 = 0) - E(9 = n/2))/N is 
of order — J c , so that the anisotropy is perpendicular. 
When the two lower sub-bands are both occupied by 
holes, the energy difference is of order — J 2 m6 2 /c because 
more holes occupy the lowest sub-band and perpendicu- 
lar anisotropy still dominates. But when holes begin to 
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FIG. 2: The minimum energy _E mm /7V and the angle versus 
p for c = 0.4 and A = 10. Inset is a plot of p versus p, revealing 
an "S" shaped curve. A Maxwell construction yields phase 
separation between the two solid circles. 



occupy the lower of the two upper sub-bands, the energy 
difference becomes positive and of order J c due to the 
larger splitting of the upper sub-bands when 9 = 7r/2. 
Assuming that J c <C Ai, the reorientation transition 
occurs close to the filling where the chemical potential 
(j, « Trp/m a b 2 crosses Ai w (7r 2 /2A 2 & 2 )(l/m; — 1/m/j). 
Very roughly, this implies that the transition from per- 
pendicular to parallel anisotropy occurs when p ~ 1/A 2 , 
independent of the Mn concentration. 

The scenario of magnetic reorientation would be re- 
versed when tensile strain overcomes the effects of carrier 
confinement so that Ai < 0. Because perpendicular or 
parallel anisotropy dominates for compressive or tensile 
strain only at very low hole concentrations, the behavior 
of the anisotropy at higher hole concentrations has led to 

the often- heard statement [ni mm that compressive or 

tensile strain is associated with parallel or perpendicular 
anisotropy. 

We now examine more closely the details of the spin 
reorientation. With both wavefunctions ipi{ z ) an( i ^(z) 
included, there are 8 hole bands rather than 4 for every 
k^ point. Since ^(CJ) = 0, the Mn spins only couple 
to the holes in ipi(z), with projection |-0i(O)| 2 = 2/L. 
The wavefunctions ipi(z) and ip%(z) are coupled by the 
off-diagonal terms in the KL Hamiltonian with matrix 



elements proportional to (n\(k x ± ik y )k z \m), which van- 
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FIG. 3: (Color online) The total energy gain of a phase sep- 
arated mixture (blue dot-dashed curve), the Coulomb cost 
associated with the domain wall (red dashed curve), and the 
total energy when both contributions are added (black solid 
curve) versus the size of the region for c = 0.5 and A = 10. 
The higher and lower values of 2Xf are associated with per- 
pendicular and parallel Mn moments, respectively. The re- 
sults below 2Xf have no physical meaning since Di and D2 
would have to be negative. 



ishes for n = m but is given by —(8i/3L)(k x ± ik y ) for 
n = 1 and m — 2. The energy E(J c ,0) of the KL plus 
exchange (KLE) model is obtained by first diagonalizing 
the Hamiltonian written as an 8 by 8 matrix in j = 3/2 
and n, m = 1,2 space. The resulting eigenvalues are in- 
tegrated over up to the Fermi level [8(. 

Fixing c and p, the energy E is then minimized with 
respect to 9. For c = 0.4, E min /N is plotted versus p in 
Fig. 2. As shown, the angle 9 corresponding to the mini- 
mum energy smoothly increases from at p\ = 0.0128 to 
7r/2 at P2 — 0.019. Of course, the chemical potential /1 
satisfies the condition dE min /dp = Nfi. We have plotted 
p versus \x in the inset to Fig. 2. This "S" shaped curve 
is typical of phase separation, where regions of different 
{9,p} coexist. The phase with high or low hole density 
is metastable so long as d 2 E min /dp 2 > or dp/dfi > 0, 
while the phase is unstable when dp/dfi < 0. Perform- 
ing a Maxwell construction for the data in Fig. 2, we find 
that phase separation occurs between regions with hole 
concentration pi = 0.0113 < p% and 9 = and regions 
with p u = 0.02 > P2 and 9 = ir/2. The phase-separated 
region is bracketed by the solid circles in Fig. 2. As the 
hole concentration increases from pi to p u , the total area 
of the regions with parallel or perpendicular anisotropy 
increases or decreases, respectively. A phase-separated 
phase is stable only if the KLE Hamiltonian is expanded 
in more than one ip n (z). 

However, phase separation of the quantum well into 



regions with low and high hole concentrations costs 
Coulomb energy, which was not taken into account by the 
KLE model. The size of the phase-separated regions will 
be determined by a balance between the cost in Coulomb 
energy and the energy gained by phase separation. The 
Coulomb and phase-separation energies are estimated by 
supposing that a charge-density wave in the Mn plane 
oscillates between fillings pi and p u within a rectangu- 
lar region in the z — plane of length D and width b. 
As described in the inset to Fig. 3, a region of length D\ 
and filling pi and a region of length D2 and filling p u 
are separated by an interface of length \p- The total 
length of the rectangular region is D = D\ + D2 + 2Xp, 
which is measured in units of the in-plane lattice con- 
stant b. For overall filling p, charge conservation requires 
that pD = piDi +p u D 2 + \f(pi +Pu)- 

To calculate the energy gained by creating a phase- 
separated mixture, we subtract the energy of a uniform 
phase with p = (pi + p u )/2, yielding the blue dot-dashed 
curve in Fig. 3 [13[. The dielectric constant for GaAs is 
used to evaluate the Coulomb energy cost, given by the 
red dashed curve in Fig. 3. The Fermi wavelength \p = 
2-K/kp is evaluated for this uniform filling p with both 
Mn orientations. While the perpendicular ^ is a bit 
smaller than the parallel fcp, the calculated total energies 
are both positive. We conclude that phase separation is 
prevented by the cost in Coulomb energy for c = 0.5 and 
that the Mn moments will be canted for fillings between 
Pi and p 2 - This is also the case for smaller Mn dopings 
which are slightly more unfavorable to phase separation. 
In the homogeneous phase with filling p, the canting angle 
9 is not affected by the Coulomb energy and 9{p) can be 
taken directly from our results. 

Since the long-range electric field contribution to the 
Coulomb energy may be suppressed in a double quantum- 
well, where the excess electronic charge on one quan- 
tum well is offset by the deficit on the other, a double 
quantum-well may exhibit phase separation rather than 
canting. By tuning the distance between the individual 
quantum wells, it may be possible to control the trans- 
formation from canting to phase separation. The depen- 
dence of the hole density on the magnetization angle will 
induce a coupling between spin- and charge-density exci- 
tations analogous to the ones observed in electronic sys- 
tems [3l ■ The emergence of phase separation will soften 
the transverse charge-density mode in a double quantum 
well. 

In Fig. 4, we plot the phase diagram of the magnetic 
orientation in the quantum well, leaving open the pos- 
sibility of phase separation. The canted region is just a 
bit narrower than the phase-separated region shown in 
the figure. The lower bound to the canted region is given 
fairly accurately by p\ « 0.012 for all c. By contrast, 
the upper bound P2 increases from 0.0135 at c = 0.15 to 
0.018 at c = 0.5. Hence, both the phase-separated and 
canted regions grow as the Mn concentration increases. 
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FIG. 4: (color online) The phase diagram of the magnetic ori- 
entation in a quantum well with A = 10, for Mn concentration 
c versus holes per cation p, showing regions of perpendicular 
and parallel anisotropy, as well as a possible phase-separated 
region. 



Due to the large size of the gap energy Ai, we do not 
expect our results to change very much with increasing 
temperature so long as T <c Ai « 29/ A 2 eV. When 
the Mn concentration c(z) is distributed along the width 
of the quantum well, our results with effective Mn con- 
centration c c s — (L/2) J _ L / 2 dz\tpi(z)\ 2 c(z) will be un- 

changed provided that (L/2) J_ L / 2 dz\ip2(z)\ 2 c(z) -C c c ff. 
We have also examined the consequences of moving the 
Mn plane from the center of the quantum well to z' = 
±L/6. Since |^i(z')| 2 = \tp2{z')\ 2 , the Mn moments in 
this plane couple equally to the holes of both wavefunc- 
tions. For c = 0.4, the canted region between p\ = 0.048 
and P2 — 0.064 is shifted substantially upwards by the 
displacement of the Mn plane. 

Although the spin reorientation transition has not been 
studied in a quantum well, several experiments have been 
performed on thin films. By measuring the remanent 
magnetization along different field directions, Sawicki et 
al. |5[ estimated that the change from perpendicular 
to parallel magnetic anisotropy in 400 nm films with 
c = 0.03 or 0.05 occurs when p w 10 20 cm -3 or 0.0045 
holes per cation. Using angle-dependent x-ray magnetic 
circular dichroism to study 50 nm films with c = 0.02 
or 0.08, Edmonds et al. 15j observed this reorientation 



transition at a hole concentration of about 5 times higher 
or 0.0225 holes per cation. 

The hole concentration in a quantum well may be con- 
trolled using a field-effect transistor, such as the one built 



by Ohno et al. [16]. Of course, the canting or phase sep- 
aration of the Mn moments would be easier to observe 
if the spin-reorientation transition happened at higher 
values of p. As discussed above, the easiest way to en- 
hance the hole filling of the spin reorientation in a nar- 
row quantum well is to displace the Mn plane by 1/6 of 
the quantum-well width. The relevant ingredients of the 
spin-reorientation transition are the single-particle gap 
Ai , the spin-orbit coupling in the valence band, and the 
magnetic coupling of the Mn ions with holes. Accord- 
ingly, we believe that our predictions extend well beyond 
the simplest confinement potential discussed here to a 
broad class of quantum-well potentials. 
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+E(p,)D 1 + E(p u )D 2 - E{p)D. 
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In the unphysical limit D — > 0, Di and D-z — > —Xf and 
Ai? PS (0) > 0. Although the intercept A£ PS (0) is propor- 
tional to Xf, it hardly differs for the two possible values 
of Xf for parallel or perpendicular moments so their av- 
erage is used to plot the blue dot-dashed curve in Fig. 3. 
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